Logarithmic current fluctuations in non-equilibrium quantum spin chains 
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We study zero-temperature quantum spin chains which are characterized by a non-vanishing 
current. For the XX model starting from the initial state | • • • TTTIII ■ " " ) we derive an exact 
expression for the variance of the total spin current. We show that asymptotically the variance 
exhibits an anomalously slow logarithmic growth; we also extract the sub-leading constant term. 
We then argue that the logarithmic growth remains valid for the XX Z model in the critical region. 
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I. INTRODUCTION 

A distinguishing feature of non-equilibrium states is 
the presence of currents 0, 0, 0| . Fluctuations of currents 
often exhibit universal behavior Q and shed light on the 
nature of non-equilibrium systems. Current fluctuations 
in classical systems have been extensively investigated, 
see e.g. 0, S i, 0, i ©, El ED, M HE US and a re- 
view [151 ] . The total current grows linearly with time and 
current fluctuations usually exhibit an algebraic growth. 
Quantum fluctuations in general, and spin current fluc- 
tuations in particular, are much less understood; even in 
the simplest systems quantum fluctuations often behave 
very differently from standard statistical fluctuations (sec 
e.g. 0.). 

In this paper we study current fluctuations in quan- 
tum spin chains. How to impose currents in spin chains? 
Perhaps the simplest way is to start with a spin chain in 
the following inhomogeneous product state 



TTTIII 



(i) 



This choice [T7| allows one to avoid complications and 
arbitrariness of coupling the chain to spin reservoirs. 

State |T]) evolves according to the Heisenberg equations 
of motion. The average magnetization profile has been 
computed for the simplest quantum chains, e.g. for the 
XX model where the perturbed region was found to grow 
ballistically fl7| . Numerical works suggest that the 
growth is also ballistic for the XX Z chain in the critical 
region (described by Hamiltonian ([2]) with |A| < 1). 

The main goal of this paper is to study the fluctua- 
tions of spin current in quantum chains. After a brief 
description of the model in section |TT1 in section IIIII wc 
present a simple derivation of the asymptotic properties 
of non-equilibrium states in free fermion systems with 
special initial conditions. In section IIVI we probe fluctu- 
ations of the current specifically in the XX spin chain. 
First we present a back-of-the-envelope calculation for 
the variance of the time integrated current; then we es- 
tablish an exact result from which we extract the long 
time asymptotical behavior. We discuss the more gen- 



eral XX Z chain in section [V] A summary of our results 
and relation to other work is given in section PVTl 



II. MODEL 



Most generally, we consider the quantum XX Z 
Heisenberg spin chain with Hamiltonian 
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Here we set the coupling constants to unity in the x and 
y directions. The coupling constant A in the z direction 
is called the anisotropy parameter. The z component of 
the total magnetization M z = X)^L-oo s n i s a conserved 
quantity in this model, and our aim is to study the corre- 
sponding current. In the following we shall focus on the 
time evolution of this spin chain starting from an inho- 
mogeneous initial state, whereby the left and right halves 
of the infinite chain are set to different quantum states 
and are joined at time zero. Such initial states provide a 
particularly convenient framework to study currents and 
their fluctuations in quantum spin chains. 

We mainly consider the special case of A = 0, where 
the model reduces to free fermions. In this system, known 
as the XX model, the time evolution can be written in a 
compact form, which enables us to perform exact calcu- 
lations. In particular, it is possible to evaluate the scal- 
ing limit of the magnetization profile and other physical 
quantities [ij], Oil • Corrections to this scaling behavior 
were considered in [H HH . 

Interesting non-equilibrium behavior was found in dis- 
ordered spin chains [12, HH and chains at finite temper- 
atures [H HH, [H| . An alternative method has been pro- 
posed to generate stationary currents in spin chains us- 
ing a Lagrange multiplier [H, [U US HI, H3| • Using this 
method fluctuations of basic quantities have been studied 
in non-equilibrium steady states in (3lj . 

The relaxation from a large class of initial conditions 
was considered in numerous studies starting from late 
sixties; see [33, K33l uyl u39| and [36fl for a review of more 
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recent work. Our focus, however, is on non-equilibrium 
states with a non- vanishing current. 



III. HYDRODYNAMIC DESCRIPTION 

We begin by describing a simple method that allows 
one to obtain the long time asymptotic behavior of a free 
fermion system by employing a continuous hydrodynamic 
description. This helps to avoid a lengthy exact calcu- 
lation and yet the final results are asymptotically exact. 
Specifically, a justification of this a ppr oach for the XX 
model is given by exact results [l7l Il9j|. which yield the 
same asymptotical behavior as the hydrodynamical de- 
scription discussed below. 

The Hamiltonian of the free fermion system can be 
written in the form 

% = Y,<^Hvk (3) 

k 

where r) k and rjk are creation and annihilation operators 
of fermions with momentum k and e(fc) is the energy of 
an excitation with wave number k. 

In the simplest situation, the system is initially divided 
into two half infinite chains, each of them being in a ho- 
mogeneous pure state. In this case, the elementary ex- 
citations can be considered initially homogeneously dis- 
tributed in each half chains. At time zero, each mode 
starts moving with velocity v{k) = e'(fc). As the ex- 
citations are entirely independent, they do not interact 
and keep moving with their initial velocities. This argu- 
ment suggests that whether an excitations is present at a 
space-time point (n, t) depends only on the ratio x = n/t. 
Moreover, keeping x fixed, a finite neighborhood of site 
n becomes asymptotically homogeneous for t — > oo. This 
physical picture is not exact due to the finite lattice spac- 
ing. However, we believe that the above description be- 
comes asymptotically exact for any free fermion system 
in the scaling limit: n — > oo,i — > oo, n/t = const. (See 
[25I ] for a rigorous derivation of this scaling limit for the 
XX model.) 

Whether an excitation is present at position n at time 
t can be decided by noting that for n > 0, the modes 
which are present were initially on the left side of the 
chain with v(k) > n/t and on the right side of the chain 
with v(k) < n/t. Similar argument applies for n < 0. 

This method can be extended to the general case when 
the two half-infinite chains are initially in mixed states, 
e.g., one can consider the situation when the two half- 
infinite chains are set to different temperatures (25|. 

As an illustration, let us calculate the magnetization 
profile in the XX model with the simplest initial condi- 
tion |T|). The spectrum of the model is e(k) = — cos(fc), 
that is, v(k) — sin(fc). Initially, all the modes are filled 
on the left, while the right side of the chain is in the vac- 
uum state. At time t, around site n > 0, the modes 
with sin(fc) > n/t are filled, that is, all the modes 




x = n/t 

FIG. 1: (Color online) Hydrodynamic description for the XX 
chain started from the | ■ • ■ TTTI-L-L " * ) initial condition. We 
consider the scaling limit where t — > 00 and x = n/t = const. 
The shaded region shows the elementary excitations that are 
present in the scaling points indexed by x. 

with ko < k < 7T — fco, where fco = arcsin(n/£). Simi- 
larly, for n < 0, only the modes with sin(fc) > n/t are 
filled, that is, the modes with — 7r < k < —n — ko, and 
k < k < 7T. For an illustration see Fig. [T] As each mode 
carries a unit magnetization, the average z magnetiza- 
tion can be obtained by simply integrating through the 
filled modes m(x — n/t) = 1/2 + j dk/2ir. This results 
in the well known profile m(x = n/t) = — ^ arcsin(:r) for 
— 1 < x < 1, and the magnetization keeps its initial val- 
ues outside this region. This limiting profile was obtained 
in [ItJ by exact calculation. 

Other applications are given in Appendix lAl 

IV. FLUCTUATIONS OF THE CURRENT 

The local magnetization current operator for a quan- 
tum spin chain can be obtained thr oug h a continuity 
equation for the local magnetization |27j. For the XX 
model this gives 

jn = s n s n+l ~ s n s n+l W 

for the current between spin n and n + 1 . (We measure 
time in units of h). The time integrated current Co, i.e., 
the net transported magnetization up to time t through 
the bond between spin and spin 1, is a quantity which 
is less obvious to define for a quantum system in general. 
However, in the case of the setup ([T]), the integrated cur- 
rent Co can be expressed in a simple way 

C , o = E« + 1 /2). (5) 

n>l 

The average of the integrated current (Co) through the 
central bond grows asymptotically as 7r _1 t [T3|; alterna- 
tively, this can be seen from the hydrodynamic picture 
described above. 
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The hydro-dynamic approach of Sec. IIIII does not al- 
low to probe fluctuations. Therefore we must return to 
the microscopic description. Below we shall focus on the 
variance of the total current D(t) = (C 2 ) — (Co) 2 . 

Let us define the left, right and total magnetization as 
follows: 



Ml = E <> 



Mr = 



n<0 n>l 

M = M L + Mr. 



(6) 



The variance of the integrated current is equal to the 
variance of the left (right) magnetization: 

D(t) = (M 2 R ) t - (Mr) 2 = (M 2 L ) t - (M L )l (7) 

Since the total magnetization is conserved and the initial 
state is an M eigenstate, the fluctuation of M remains 
zero for any time t 

{(M L + M R f) t - (M L + Mr} 2 = 0. (8) 

By exploiting this property we can rewrite J7]) as 



In this expression, the factor 1/2 appears since correla- 
tions between evenly spaced sites vanish. 

This argument remains valid for a more general class 
of initial conditions, where on the left (right) half of the 
chain the fermions are filled up to the Fermi energy fi^ 
(hr). In this case, the asymptotic state — which builds 
up near the origin — includes fermions with momenta 
varying from — fcjj to Icl, where fci and kR are the Fermi 
momenta corresponding to [Xl and fiR. (For an illus- 
tration see Fig. [3J more details are given in Appendix 
IA 11 ) The correlation function p z (n) for these asymp- 
totic states can easily be calculated [23]. One finds that 
in general it behaves as p z (n) = — a 2 sin 2 (nip) , where 

= (fcL + fcfl)/2 [27j. As the <p dependence of the asymp- 
totic form averages out we conclude that the result (fT3|) 
is unchanged for this class of initial states. 

The exact evaluation of D(t) is based on ([9]). Follow- 
ing the strategy of [13] we write sf in terms of the local 
fcrmionic creation and annihilation operators c 1 , c as 



c 



1 



(14) 



D(t) = (M L )t(M R ) t - (M L M R ) t 

= E « s ?>*« 

K0,ro>l 



z .s z m )t). (9) 



In the Heisenberg picture, the time dependence of these 
operators, under the dynamics of the XX chain, has a 
simple form 



Before presenting an exact calculation we provide a 
back-of-the-envelope derivation of our main result. The 
idea is to evaluate D(t) by substituting correlations in © 
with their stationary values in the local state which builds 
up at the origin for t — > 00 (see [17J and our Fig. [J) . The 
reason is that the main contribution comes from those 
spins for which (I — m) is not too large, and for t ^> 1 
these "points" are located near the origin. This leads to 



D(t) = -J2np z (n), 



n>0 



where 



p z (n) 



'fc^fc+n/ 



{sl)(s 



k+nl 



(10) 



(11) 



In the (homogeneous) "maximal current" stationary 
state of the XX model, the correlator p z (n) takes the 
same form as in the ground state [27j where it is given 
by the well-known expression: 



c n (t) 



E 



(15) 



where J n (t) are the Bessel functions. Inserting this into 
one gets 

D(t)= e E *- a+p - i+s 

i<0,m>l a,/3,7,<5 

x Ja-l(t)Jf3-l(t)Jy-m(t)Js-7n(t) 

X ((4 C /3) ( C 7 C <5> _ ( C 1 C P C \ C S)) ■ (16) 

The expectation value in the above formula is taken in 
the initial state. One finds 



-S a ,s<>i3 n Hat, fil 







otherwise 



(17) 



p z (n) 



1 







n 
n 



odd 
even 



(12) 



This gives a logarithmical divergence for D(t), which one 
can regularize by truncating the sum in (|10[) . For a finite 
but large t the volume of the region around the origin, 
which can be described by this maximal current state, 
grows linearly with t. Hence we choose the upper limit 
in the sum in (|10[) to be proportional to t and obtain 



2tt 2 



ln(t). 



(13) 



for initial states which are product states of individual 
spins pointing either up or down (like (H}). Here, a |, 
f3 [ are shorthand notations for a, (3 with s a =1,sp =j. 
Using identities J- k (t) = (-l) k J k (t) and [46| 



E Jk+p(t)Jk+q(t) 



fe>l 



J p (t)J q+1 (t) - J p+1 (t)J q (t) 

2(p-q) 



for the sums over I and to, one obtains 



at, PI 



J a -l(t)Jp(t) - J a (t)Jfi-l(t) 



a-0 



(18) 
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This expression for our initial condition (JXJ) becomes 



XXZ SPIN CHAIN 



D(t) 



4 



E 

l,m>l 



Jl-l(t)J m -l(t) + Jl(t)J m {t) 

l + m-1 



(19) 



This is an exact expression for the variance of the current 
which is valid at any time t > 0. From this formula we 
have deduced the long time asymptotic behavior 



D(t) = -L(lnt+C9 



(20) 



with constant C = 2.963510026 . . .. The leading term co- 
incides with our heuristic argument p3[) . The derivation 
of (|2TJ|) is relegated to appendix iBl 

We also evaluated expression (fT9|) numerically and 
plotted it in Fig. [5] for t < 50. We observe a logarith- 
mic increase in time. The numerical estimate for the 
constant, C w 2.9633, is in a good agreement with the 
exact result. 

In addition to this logarithmic growth, one observes 
oscillations with decreasing amplitude. We found that 
the formula 



D(t) 



1 

2^2 



hit 



(21) 



with C w 1.95 (and C given exactly in (|B23[) ) gives a 
very good fit to the numerical data even for relatively 
short times. Since the mean current is 1/tt, in aver- 
age one fermion crosses the origin in time it. Hence the 
cos2i oscillations can be interpreted as a consequence 
of the quantum nature of the magnetization: each pass- 
ing fermion causes a bump in the fluctuations. Similar 
arguments were used to explain oscillations in the mag- 
netization profile in the same system [2~H . 




t 



FIG. 2: (Color online) Shown is the variance D(t) vs. time 
t. Red curve is a result of the exact numerical evaluation 
of (|19[) . blue curve shows the result p0[) of our asymptotical 
analysis. On top of the logarithmic growth we find a sub- 
leading oscillating term. Eq. (|21[) gives a fit, which is almost 
indistinguishable from the numerical data (red line). 



In a general XXZ chain ([2]), the coupling is nonzero 
in the z direction. There is no explicit solution for the 
time evolution in this general case. Due to the symmetry 
in the XY plane, however, the z component M z of the 
total magnetization is still conserved, hence the magne- 
tization current can be studied. It has been investigated 
numerically by Gobert et. al. [l|| . 

Our main interest lies in the so-called critical region 
(—1 < A < 1); we shall also discuss the isotropic ferro- 
magnetic (A = 1) and anti-ferromagnetic (A = — 1) spin 
chains. The XX model belongs to the critical region, and 
the behavior in the entire critical region is believed to be 
similar to the behavior of the XX model. In particular, 
the magnetization profile plausibly scales linearly with 
time: m(n,t) — ► M(n/t) [lq . On the other hand, due 
to the existence of kink- like ground states [33, [38| , the 
magnetization profile is expected to become frozen when 
|A| > 1 [HI]. Algebraic scaling m(n,t) — > M(n/t a ) seems 
to emerge for the XXX model (|A| = 1) in the scaling 
regime n — > oo, t — > oo with n/t a kept finite. Numeri- 
cally the exponent is a = 0.6 ±0.1 [18|, so the non-trivial 
part of the profile is sub-ballistic. 

For the XX spin chain, we have obtained the correct 
current fluctuations in the leading order from the simple 
formula (JTUJ) when the upper limit of the sum was chosen 
to grow linearly with time. The reason for this choice 
of upper limit is that the front and the whole profile 
"moves" linearly with time, hence the cutoff must behave 
similarly. We shall use (fTU)) also for the XXZ chain. 
We shall assume that the upper bound moves linearly, 
namely as vt, in the critical region (—1 < A < 1). The 
actual value of the 'velocity' v is unknown, but it does 
not affect the leading order term anyway. 

The next issue is whether one can use the equilibrium 
spin correlations p z {n) in the presence of current. For 
the XX chain we know [13, H3] that current does not 
affect the z component of the correlations significantly 
(may introduce a modulation), only the x and y compo- 
nents. Here we boldly assume the same for the XXZ 
model, at least for its large distance behavior. Thus we 
use the equilibrium correlations in (|10p . The asymptotic 
formulae for p z (n) are [39j, [40j 



p*(n) 



—8n 2 
-[l + (-l)"](2^ 2 )- 1 n- 
{~l) n An-^ 2s ~5n- 2 



< A < 1 
A = 
-1< A < 



(22) 



A 



-1 



where we used the shorthand notation 

1 



(23) 



The amplitude 
A = A(A) has been guessed relatively recently (see [39(), 
yet we do not need this result. Indeed, the leading oscil- 
lating term in p z {n) has the exponent a — 4ir 2 5 varying 



47r arccos(A) 
We now insert (|22|) into equation (jXQ 
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in the range 1 < a < 2 when the anisotropy parameter 
varies in the — 1 < A < range. Because this leading 
term oscillates, we form pairs and find that the oscillating 
terms in (jTUJ) yield the contribution 

An- a+1 ~ A(n+l)- a+1 ~aAn- a (24) 

that decays faster than n _1 . Hence the oscillating term 
provides merely a constant contribution to the variance 
D(t) while the sub-leading n~ 2 term results in the leading 
logarithmically diverging contribution. Therefore 

D(t)=6hxt (25) 

This prediction implies that the logarithmic behavior of 
the variance is universal. The amplitude diverges at 
A = 1; analytically S -> (4tt)- 1 [2(1 - A)]" 1 / 2 as A t 1. 
This divergence is not very surprising since the isotropic 
Heisenberg ferromagnet apparently exhibits a truly dif- 
ferent behavior. 

In the other extreme A J, — 1, the amplitude S ap- 
proaches the finite value (47r 2 ) -1 . However, as indicated 
by the last formula in l2"2"j) . the oscillating asymptotic 
is Bn^ 1 \/\an for the isotropic anti-ferromagnet [4lj . A 
calculation similar to ([24)) gives n _2 \/lnn after canceling 
the oscillations. This leads to the following quite surpris- 
ing behavior: 

D~ [ — VIn^~ (lni) 3 / 2 (26) 
J n 

Note that for the isotropic anti-ferromagnet the magne- 
tization is non-trivial in the interval that grows slower 
than linearly with time. (The natural guess is the diffu- 
sive \ft growth.) However, the upper limit in the integral 
in (|26p would affect only the pre-factor. 

Thus our tentative predictions for the variance of the 
time integrated current are: (i) The enhanced logarith- 
mic growth D ~ (hit) 3 / 2 for A = —1; (ii) The universal 
logarithmic behavior D (t) = 5 hit with known pre-factor 
(|2"5)) in the critical region — 1 < A < 1. 

VI. DISCUSSION 

We have studied the fluctuations of the time integrated 
magnetization current in the quantum XX chain that 
evolves starting from an inhomogeneous non-stationary 
initial state fl}. We have derived an exact formula (fl9|) 
for the variance D(t) of the current. We have shown 
that the variance increases logarithmically in the long 
time limit (|20[) . which is consistent [42J with numerical 
evaluation of the exact formula (see Fig. [2]) . In addition 
to this logarithmic growth, we have observed oscillations 
with decreasing amplitude. We have argued that this 
logarithmic leading order behavior remains unchanged 
for a more general class of initial conditions (where the 
magnetization on the two half-chains is not saturated). 

These small logarithmic fluctuations reflect the ideal 
conductor nature of the integrable XX quantum chain. 



The current simply "slides" through the system ballis- 
tically with no disturbance, hence the tiny fluctuations. 
Conversely, if an impurity is present at the origin, the 
variance grows linearly with time J43|. Similarly, in 
stochastic particle systems 0, H, @, 0TB [1] the noise - 
which is intrinsically present in these models - generates 
algebraic fluctuations. 

We have argued that current fluctuations in the in- 
homogeneous XXZ model are also logarithmic (in the 
critical region) . Our arguments are heuristic and a more 
rigorous derivation is a key challenge for future work. In- 
triguingly, fluctuations seem more tractable than e.g. the 
average magnetization profile in the XXZ chain, which 
is completely unknown. 

For free fcrmion systems, current fluctuations were 
found to be asymptotically Gaussian [43[ and therefore 
the variance provides a complete characterization of the 
full current statistics. Moreover, according to [44| this in- 
dicates that the entanglement (between the left and right 
halves) is simply proportional to the variance of the cur- 
rent (with a factor 3/7r 2 ). In order to check this relation 
we compared our results for D(t) to numerical results for 
the time dependent entanglement entropy (for the same 
model and initial condition) presented in Fig. 19 of [l8[, 
and we found a good agreement in the leading order. 

Little is known about higher moments of current fluc- 
tuations for interacting fermions. There is no reason to 
believe that they are Gaussian and the full current statis- 
tics is very difficult to probe (both theoretically and ex- 
perimentally) for interesting interacting fermion systems 
|45j |. (Even for classical interacting particles, the deriva- 
tion of the full counting statistics is usually a formidable 
challenge, see e.g. @, 1 1 H H [H, [l|.) For the 
XXZ model, however, it could be possible to compute 
higher order cumulants by employing the heuristic ap- 
proach which we have applied to computing the variance. 
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APPENDIX A: HYDRODYNAMIC 
DESCRIPTION: EXAMPLES 

1. XX model 

In section Mil we demonstrate how the hydrodynamic 
approach works for the XX model in the simplest case, 
i.e., when all the spins point up (down) on the left (right) 
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side at t = 0. One can easily recover the exact results 
for a slightly more complicated initial state as well. Con- 
sider a product state composed of a ground state (in an 
appropriate external field) with ttll {itir) magnetization 
on the left (right) side of the chain. That is, the modes 
-kL(R) < k < k L(R) with k L ( R ) = vr(l/2 + m L(R) ) are 
filled on the left (right) side. Without any restriction, in 
the following we assume that mj, > m R . 

There are two qualitatively different cases. (1) If the 
sign of rriL and m R are equal the scaling profile consists 
of three segments. Here we assume that tul > 0, m R > 
but the case mj, < 0, m R < is entirely analogous. Ap- 
plying the argument of section HTTl hcre the magnetization 
profile takes the following scaling form: 



m L +m R 

mix) = { ** + \ 
m R 



arcsin(x) 

27T • 



< x < v L 
v L < x < v R 
v R < x, 



(Al) 



v(k R ) 



with vl = w(fci) = cos(7rmi) and v R 
cos(7rmfl). and m(—x) = rriL + m R — m(x). 

(2) A qualitatively different case is when tol > and 
m R < 0. The argument of section HTT1 leads to the follow- 
ing asymptotic profile: 



m(x) 



m L +m n 
2 

mi? i 1 arcsinfa;) ^ - 

-f" + I 2F^' v L <x< VR 

mR + l - HSiilM 

m R 



< x < v L 



V R < X < 1 
V R < X, 



(A2) 



As an illustration see figure[3] Here we assumed vl < v R , 
as the opposite case is entirely similar. In all cases m(x) 
has the symmetry property m(—x) = rn^ + m R — m{x). 
We mention that the special case of m R = —rriL was 
considered in detail in [l?], Hfl . 



2. Dimerized XX model 

In all the above example the speed of the front was 
one, but it is not necessary. Consider, as an example, 
the dimerized XX model [3a| 




-1 -v R -v L v L v R 1 
V = n/t 

FIG. 3: (Color online) Hydrodynamic description for the XX 
chain started from an initial condition where the left (right) 
side of the chain is in the ground-state with magnetization 
rriL (mR). We consider the scaling limit where t — > oo and 
y = n/t = const. The shaded region shows the elementary 
excitations that are present in the scaling points indexed by 

y- 

5. That is, the front still moves ballisticly, but it can go 
arbitrarily slowly. By inverting the function ui'(k) = x/2 



ki 2 = arccos 



-x 2 ± y/(l~x 2 )(S 2 



(A5) 



with < x < 6, after some elementary calculations, one 
obtains the magnetization profile as well 



m(x) = 1 (&2 — hi) = arcsin 

2 27T 7T 



Surprisingly, the dimerization only rescaled the asymp- 
totic magnetization profile, and the result agrees with a 
homogeneous chain with coupling S at each bond. That 
means that the weaker links behave as bottlenecks and 
they govern the dynamics of the chain. The average time- 
integrated current through the origin is tS/n, which was 
also found to be linear in time in [1 81 ] . The value of their 
current agrees with our result (see figure 2]) . 



H — Jm(s m S m+1 + Sm s m+l) (A3) 

with J m — 1 for odd m, and J m = S < 1 for even m. The 
spectrum of the model is 

uo(k) = ±-yT + <5 2 + 25cosfe, (A4) 

where k is a sublattice wave number going from — n to 
7r (an excitation with sublattice wave number k has mo- 
mentum k = k/2). The two signs in (|A4[) indicate two 
branches with k 6 (— 7r/2, 7t/2). For the simplest ini- 
tial condition with mj = — m R — 1/2, the fastest mode 
present on the left [at k = ±(ir — arccos 5)] has speed 



APPENDIX B: LARGE TIME ASYMPTOTIC OF 

D(t) 

Here we study the large time behavior of (|19p. As 
the summation in (|19|) goes to infinity, we use a general 
large t asymptotic formula [461 ] of Bessel functions which 
is valid for all m < t: 

fY cos tf(ji) - j /jx 

J "'" = V;n a- ».)■/' +0 (i)- < B1) 

where 

/(M) = V 7 ! - M 2 - M cos" 1 fj, , M=-r- ( B2 ) 
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FIG. 4: (Color online) Shown is the average transferred to- 
tal magnetization M + = (52 n>0 (S„ + 1/2)) as a function of 
time started from the | ■ ■ ■ T TT 4-4-4- ' ' ' ) initial condition in the 
dimerized XX model. The data points show the result of 
a time-dependent DMRG calculation published in [lg| . solid 
lines indicate the prediction of the hydrodynamic description. 
The coupling constants for the weak link are 1, 0.8, 0.6, 0.4, 
0.2 (from above to below). 



We will also use the scaled variable A = l/t. For m > t 
the Bessel function J m (t) becomes exponentially small 
thus we can safely neglect those terms in (|T5| . Hence in 
the double sum we consider only the range A,/i < 1. 

Based on (IB1[) the asymptotic form of a Bessel func- 
tion with a shifted index can be obtained by simply Tay- 
lor expanding /(//) in the argument of the cosine (the 
term coming from the expansion of the denominator is 
negligible), which gives 



7 M AT^cosa- \]\- sin a 

J m -x{t) = y {1 _ fi 2 ) i / 4 + 



7 (B3) 



where we introduce the shorthand notation 

7T 7T 

4 ' fc = */0*)-4- 



a 



Now using the asymptotic formulas (|B 1 1) and (|B3[) in 
T9|) we obtain 

{2/irt)T 



Jl-l(t)Jm-l(t) + Jl{t)J m {t) 



where 



(1- A 2 ) 1 /4(l- M 2)l/4' 

(B4) 



T = (1 + A/i) cos a cos 6 + \J\ — A 2 \/l — /i 2 sin a sin 6 
— A\/l-/j 2 cos a sin b — — A 2 sin a cos fe. (B5) 

From hereafter we omit noting the relative 0(1 /yt) cor- 
rections. 

When \,n < 1, we can replace (1 - A 2 ) 1 / 4 (l - /J 2 ) 1/4 
by one in (|B4[) . and the expression for T simplifies to 



T = cos(a - b) = cos[t/(A) - tf(fi)]. 



Using (|B2|) we expand /(A) and /(//) to give 



tf(X)-tf(u-X) = 



tt('-to) I 2 



2t 



When I, m <C y/i, we can keep only the first leading term 
of this sum. In this limit (IB4I) reads 



Jl-l(t)J m -i(t) + Ji(t)J m (t) = — cos 

7Tt 



7T(Z 



(B6) 



Now in (HHJ) we write X^,m>i 
divide the summation through n into four parts 



En>2 Ej-|-m=Ti! ano - 



[t l/2-« 

E = E 

n>2 n=2 



[*V»+.] [t ] 

E + E 

n=[tVa-*]+l ra=[tV2+']+l 



E ( B7 ) 

n=[t] + l 



D(t) = Di(t) + D 2 (t) + D 3 (t) + D 4 (t) 



where [•] denotes the integer part. The corresponding 
terms in (Tl9t will be referred to as 



(B8) 



respectively. We choose the exponent in (|B7[) in the range 
< e < 1/2; later we shall take the e — > limit. Since 
ji/2-e ^- ag ^ _^ ^ m £) 1 (^) we can use the simple 

formula (|B6[) . and immediately perform one summation. 
Thus the first sum becomes 



[t 1 



Dx(t) 



E 



i 



n -V 



which in the t — > oo limit is simply 

(1/2 - e) ln(i) + 7B + ln(2) 



Di{t) 



2tt 2 



(B9) 



(BIO) 



where -fE = 0.5772 ... is the Euler constant. 

Now we consider the contribution from the region 
f i/2-e < n < £i/2-e. Using cos 2 (-) < 1 we find that 



D 3 (t) < 



[t l/2+« 

E 



i 



2c 



n=[*i/2-«] + l 



n - 1 



ln(<) 



(Bll) 



for fixed positive e and large t. Hence one can see that for 
small e the contribution from ^(i) becomes negligible 
(as compared to D\{t). 

In D$(t) and D±{t) we can replace summation by in- 
tegration, since t 1 l 2+<L — > oo as t — ► oo. We take the 
square of (|B5j) and drop rapidly oscillating terms (like 
sin a cos a) while sin 2 and cos 2 are replaced by 1/2; this 
results in replacement of T 2 by (1 + A/i)/2. Overall, we 
find that 



Ds(t) 



1 

27 2 



[*] 

E 

n=[tV2+«]+l 



1 



(n-1) 2 



E 



1 + A/i 



l n V(l-A 2 )(l-/i 2 ) 



(B12) 
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and 



1 [2tl 1 



(n - l) 2 

1 + A^ 
„V(1-A^)(1-^)' 



ra=[t] + l 

X V ; ^ IBW) 

^ ,/(l- A 2 Vl 77 



Note that there is only an exponentially small contribu- 
tion from n > 2t terms. 

Notice that the second sum in (|B12IIB13| can be re- 
placed by 



1 



n Jo 



F(v, x)dx, 



where we used the shorthand notation 
1 + v 2 x(l - x) 



F(v, x) 



y/[l ~ V 2 X 2 ][l - V 2 {\ - x) 2 ]' 



(B14) 



(B15) 



and introduced v = n/t = A + /x and x = \jv. For 
t — > oo the first sum in (|B12IIB13)) can also be replaced 
by an integral, which leads to 

D ^) = ^l — I dxF{v,x) (B16) 
2tt z J t ,-i/2 v J 

It is useful to rewrite (|B16|) as 
dv f 1 



(B17) 



2?H J t e-l/2 V J 

, (l/2-e)lnt 
+ 2^ 2 



where the first integral is convergent in the limit t — > 
oo,e — * 0. Therefore the contribution gathered in the 
£i/2+£ <c n < t region is 



D 3 (t) = 



(l/2-e)lnt + C 3 
2tt 2 



(B18) 



where 



f d f 

C 3 = — dx[F(v,x)-l] = 0.34929294... (B19) 
Jo v Jo 



By a similar argument the contribution (|B13|) from the 
t < n < It region remains finite in the t — > oo limit: 



D 4 (t) 



C4_ 
2tt 2 



(B20) 



with 



dv r 



C 4 = / — / dxF(v,x) = 1.34385423... (B21) 
Ji v Jx_ v -i 



Combining the contributions (|B10|) . IjBlip . (|B18|) . and 
(IB20[) from the three regions of (|B7[) . and taking now 
the e — ► limit, we obtain 



(B22) 



where 

C = 7e + In 2 + C 3 + C 4 = 2.963510026 . . . (B23) 
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